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We have studied the thermodynamic properties of a multiferroic that couples ferromagnetic and 
ferroelectric order. Some of the results are independent of the form of the free energy. We calculate 
the temperature dependence of the electric, magnetic, and magnetoelectric susceptibilities. The 
cross susceptibility has a temperature dependence related to the mixed (with respect to E and 
B) derivatives of the specific heat. The phase transitions are all second order. In particular, the 
phase boundary Tm(E), where Tm is the lower magnetic transition as a function of electric field, 
is described by the Ehrenfest relation. The magnetoelectric susceptibility is nonzero only below the 
lower of the two transition temperatures. We study the properties of the specific heat, with and 
without the inclusion of gaussian fiuctuations. The perturbative renormalization group is used to 
understand the fixed points of the theory, and we include a discussion of the effect inhomogeneities 
have for this model. 

PACS numbers: 

I. INTRODUCTION 

Multiferroics couple polarization, magnetization, and the elastic response in a material. The subject has been 
considered interesting for some time now. Interest in the electric control of the magnetic properties of a system for 
appHcations requires the magnetic and polar orders be coupled as strongly as possible. There are materials that 
exhibit multiferroic properties and have been subject of several reviews recently^'^'i^. Many physical systems can be 
described by the interplay of more than one order parameter, but it is the new context of multiferroics which leads 
back to this classic problem. 

Dzyaloshinskii^ first pointed out that, on symmetry grounds, the coexistence antiferromagnetic and ferroelectric 
order in Cr203 was possible. Astrov^ confirmed the prediction for Cr203 by measuring the magnetization due to 
an electric field. Rado et aS specifically demonstrated the nature of the interaction between the magnetic and 
electric degrees of freedom by reporting that the temperature dependence of the cross susceptibility (defined below) 
was different depending on the relative orientation of the magnetic field and the staggered magnetization in the 
antiferromagnetic Cr203. During the 1960s and 1970s several multiferroic materials were discovered, such as BiFeOa, 
BiMnOa, and the boracite Ni3B20i3l. A review by Smolenskii and Chupis^° is a compilation of the subject until 
the early 1980s. Recently, the discovery of a new class of multiferroics, notably TbMnO^^, reignited interest in 
multiferroic research. A broad class of these multiferroics follows the chemical formulas, RMn03 and RMn20g^ , 
where R is a rare earth ion. The appearance of spiral magnetic order causes a polarization. In this more recent class 
of multiferroics, the ferroelectricity is referred to as improper meaning that the primary order, the magnetism, is said 
to induce the electric polarization which occurs as a consequence of the symmetry. Mostovojsi^ provided a general 
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argument based on symmetry and Ginzburg-Landau theory for why this happens. This is to be contrasted with proper 
ferroelectrics, where a spontaneous polarization exists of its own accord. We may then classify the macroscopic view 
of this behavior as follows: in an improper ferroelectric Tm and T^; coincide (P occurs as soon as the appropriate M 
does) while in a proper ferroelectric T^; is independent from (usually greater than) Tm. This paper is a study of the 
later case, its properties and fixed points of the theory. 

The physics of two coupled order parameters has been studied in various contexts, see 18123 -25 and references 



therein. More recently, the dynamics and domains in ferroelectric ferromagnets have also been studied2ii2^i2i. We 
consider the thermodynamic response of a system which does not depend on any specific form of the free energy, then 
specialize to ferroelectric-ferromagnetic order to investigate the temperature dependence of the susceptibility which, 
to our knowledge, does not appear in previous work. An outline of our paper is as follows: we begin with a catalog of 
thermodynamic results including Maxwell relations and a bound on the magnetoelectric effect^. In the next section, 
we specify the simplest form for the free energy and examine the susceptibilities at both transition temperatures. We 
find that in the absence of a fluctuation-induced response, the magnetoelectric susceptibility is only nonzero when 
both P and M are nonzero. Finally we examine the fixed points for the theory and the effect of inhomogeneities. 

II. THERMODYNAMICS 
A. Maxwell Relations 

The Helmholtz F(P,y,M,T) and the Gibbs G(E,P,B,T) free energies are functions of the volume V, magneti- 
zation M, temperature T and Electric field E, polarization P, pressure P, and magnetic field B. They are described 
by, 



dF = -SdT - PdV + EdP + BdM (1) 

dG = -SdT + VdP - PdE - MrfB (2) 

These are the definitions of thermodynamic observables as derivatives of the free energy. By noting the equality of 
cross derivatives, it follows that, 

9e'^-9b'^ (3) 
Which is the reciprocity statement for the cross susceptibility Xi2 in Eq. [5T] By looking at the other cross 
derivatives, the following relations can be derived: 

9E2 dT^ dB^ 

and 

(^) 



dBdE 9T2 



Here x is the differential susceptibility (the subscript E for electrical and M for magnetic) , k, is the compressibility, 
V is the volume and C the specific heat at constant volume. 

The mechanical stability of the ground state also leads to the inequality 

d^F d^F , d^F 



(— ) < (7) 



which in turn implies that 

X?2 < XeXm (8) 

This inequality was first derived by Brown et al^ and remains an essential testing ground against unsupportable 
approximations . 

The Ehrcnfcst definition of the order of a phase transition needs to be considered carefully where there is more than 
one mechanical field. The thermodynamic phase boundaries for a first order transition and a single order parameter 
are described by the Clausius-Clapeyron equation. For a second order phase transition, the corresponding equation 
is the Ehrenfest equation^'^. 

The Clausius-Clapeyron equations for Tm(E), which can be obtained from dGi — —SidT + VdP — PdE — MdB 
by choosing two points on either side of the critical line in the T-B plane, takes the form: 

dTu Pi - P2 



9E Si - S2 

For completeness, we present all four Ehrenfest equations with the two mechanical fields, E and B: 

dTES" TeAxe {OTmY TmAxe 



(9) 



9E y AC ' V 9E y AC ^^^^ 

dB J AC ' \dB J AC ^^^^ 

In a ferromagnet (or ferroelectric) in any finite field there is no phase transition. And yet, if the transition is second 
order at a finite field, as it can be for an antiferromagnet, the evolution of Tc should be described by the Ehrenfest 
equations above. When the two order parameters are coupled, as we will show below, Axm is zero at Te but Axe 
is non-zero at Tm- The ferromagnetic transition at rM(B,E) is second order in E, unlike the uncoupled result. 

B. Adiabatic Processes 

Adiabatic processes are used in cooling, in a Joule- Thompson decompression or demagnetization at low temper- 
atures. In a multiferroic, the entropy S(E, B, P, T) depends on the external fields. The cooling arising from an 
adiabatic process is described by dS=0 which upon substituting the Maxwell relations becomes: 

^rfE+^^dB + ^dP+^dT = (12) 
dT dT dT T ^ ' 

^ = -^^ (13) 
dE C^dT ^ ' 

We will defer any explicit calculation of the integral for the which depends on the explicit temperature dependencies 

of the specific heat and also of the order parameter. 
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III. FREE ENERGY FUNCTIONAL 

The free energy of a system with more than one vector order parameter can be written as: F — Fe + Fm + Fint 

Fe = ^ + 6bP^ (14) 

Fm = 7, + &mM^ (15) 

^Xmo 

F, = fc(M2p2) (16) 

XbJ = aso - 1^ (17) 

Xmo = - 1^ (18) 

We are considering here a ferromagnet described by Fm and a ferroelectric represented by -Fe. The interaction 
between the two degrees of freedom is given by Fi . The possible inhomogeneous instabilities will be considered in the 
summary section below. The coefficients of the free energy here are all constants with the exception of the various 
susceptibilities which are assumed to be described by the usual Curie- Weiss law for localized moments. For specificity 
we take Teq > Tm as the transition temperatures respectively for the electric polarization and the magnetization. 
In general, the interaction between the magnetization and electric polarization must be a scalar and could be in the 
form fciM^P^ + fc2(M.P)^ which is required by the time reversal invariance. This term determines the relative angle 
between M and P. Once the relative angle is determined, the free energy becomes of the form in EqlTHlwith its 
coefficients renormalized. There is another motivation for considering a biquadratic interaction of this form, beyond 
it being the simplest consistent with symmetry. If we think of the exchange for local moments and expand, since the 
dipole moment is classically proportional to the displacement, the coupled elastic and Heisenberg interaction, through 
\7J{x^j)PA'P, result in a P^M^ form. While this argument is only heuristic, it informs the choice of free energy. 

The equations of state are the conventional thermodynamic equations for a Helmholtz energy: 

dF 

E=aj; (19) 

dF 

- am P») 

For a homogeneous system, the free energy above is quite sufficient to derive all thermodynamic properties. In a 
multiferroic, the linear response includes cross susceptibilities Xi2, defined by 



5M^ XMB + X12E ^^^^ 
SP= X12B + X-EE 



The free energy functional in the conventional form includes the external fields — P.E and — M.B, and thus the 
ground state is a minimum of f. The dimensionless free energy / — F/Fe{0), is written in terms of: p — P/Pq, 

E 



m = M/Mo, I = Fm(0)/Fe(0), e = ™Po, ^-mi = B.M/F^e(0). 
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f^^ = fe + frn + ^km^p" ^ 2(-^ - l)p' + / - p.e + l{2{^ - l)m2 + - m./3) + 2k{m'p') (22) 

-r_BO J BO J MO 

The equations of state in the scaled variables becomes, 



i = {^-l)p + p^ + kw?p (23) 



^ = (-^-l)m + m3 + ^mp2 (24) 

4 J MO ^ 



A. Order Parameter 



We consider Teo > Tmo- The solutions are: 



Po = "lo = T > Teq > Tm 

ml = 0, pI{T) = {1-^) Teo>T> Tm. 

ml{T) = -^(1 ~ ^), pI{T) = -i^(l - ^) T^^o > Tm > T 

Teo l^k'J^ Tmo i-^Imo 

-'a/o -'bo 

As shown in Fig. [1] the mean field order parameter P appears at Teo continuously. 




0.0 ^ ' ' ] 
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FIG. 1: (Color Online) The polarization and magnetic moment versus temperature. P shows a jump at Tm when M condenses. 

Since we have chosen Teo/Ta/o > 1, we have Te/Teo > 1 but Tm/Tmo < 1- The electrical transition takes place 
at Teo, which is unchanged by the interaction, and Tmo is re-normalized by interactions. The physical magnetic 
transition takes place at Tm, not Tmo- As shown in Fig. [1] there is a kink in po at Tm- This kink is due to the fact 
that the scale temperature for po, for T < Tm is Te > Teo- 
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B. Susceptibility 

Following the equation of state, we can derive the electric, magnetic, and the cross susceptibilities as defined in 
Eq. [211 They both take a Curie- Weiss form, also have a structure at the other transition temperature. Figures 
(I2EI) show these features for a specific choice of the parameters(k= .4,j = .3). The cross susceptibility diverges 
at Tm- As expected, it vanishes above Tm- This excludes the possibility of a fiuctuation induced response, which 
we neglect in our study. The ferroelectric susceptibility, exhibits a divergence at its transition temperature and an 
anomaly at the magnetic transition temperature. The magnetic susceptibility demonstrates paramagnetic behavior 
above it's transition and diverges at Tm. It shows a cusp at the ferroelectric transition temperature coinciding with 
the smooth onset of P. The inverse electric susceptibility, Xe^ i is zero at its transition temperature and jumps at the 
magnetic transition temperature. The size of this jump depends on the energetics of the ferroelectric and magnetic 
free energies, and occurs because both D and the numerator vanish as T approaches Ta/ from below. The inverse 
magnetic susceptibility, Xm ' is zero at Tm , and shows a change at the ferroelectric transition temperature due to the 
onset of P, resulting in a cusp. 
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FIG. 2: (Color Online) (Left) The ferroelectric susceptibility, the blue solid line, and the magnetic susceptibility, the dashed red 
line. (Right) The inverse susceptibilities, the blue solid line, Xb^i s^nd the thin dashed line representing the inverse magnetic 
susceptibility, Xm ■ An auxiliary dashed line is also present to underscore the slope change in Xm Teo- 

The formal expressions are, 



XEiT > T, 



EO) 



4(^-1) 



, Xe{Tm <T < Teo) = 



Teo' 



(26) 



XEiT < Tm) - ^((^ - 1) + + ^pI) 



(27) 



Xm{T > Tm) 



1) 



(28) 
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FIG. 3: (Color Online) The cross-susceptibility becomes nonzero only below Tm and has the opposite sign of the other 
susceptibilities. 



with 



Xm{T < Tm) = - 1) + 3Po + J^o) (29) 

1 k 

Xi2 — (-2-moPo) (30) 
n If 1^ , q 2 , , 2ir. T ,^ , o 2 , ^ 21 (2fcmopo)^ 

^ = [(^ - 1) + 3po + fc™o] [(^ - 1) + + -^Po] -g (31) 

The results have been plotted in Figs. ([2][3]) for certain specific values of the parameters representative of the general 
behavior. These parameters include the two bare transition temperatures(T£;o=39K, T£o=28K), the interaction 
parameter k and the free energy ratio I. In all of the results discussed so far, fc^ < The ground state is different 
in the opposite case and leads to a correspondingly different phase diagram. The cross susceptibility Xi2 is non-zero 
only below Tm- 

C. Specific Heat 

The free energy describes two second order phase transitions which are coupled. In a mean field like analysis, that 
corresponds to two discontinuities at Tj\/ and T^o- The algebraic results are: 

^ = T>Teo, ^ Tm<T<Teo, T{-^){-^ + -L- ) t<Tm (32) 



IV. INHOMOGENEOUS EFFECTS 

There are three gradient terms in the free energy, one for each of the order parameter as well as one in the interaction 
term. In the following we first consider the stability of a uniform ground state with respect to inhomogeneous 
perturbations. This is followed by a discussion of the fluctuations. 
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FIG. 4: (Color Online) The specific heat with the same parameter choice as the susceptibility. 

A general interaction between P and M, involving space gradients (and a scalar with respect to inversion and time 
reversal) can be written as PiMjMkdi. The gradient can operate on either M or P. The four indices (i,j,k,l) should 
be contracted to turn the interaction in to a scalar. Thus the interaction terms are 

ziM.{{M.V)P) + Z2Af^(V.F) + Z3(P.Af)(V.M) + ZiP.[{M.V)M] (33) 

Since all terms are quadratic in M, there is a shift in "Tm proportional to the gradient of P. The first pair of terms 
can be integrated by parts and turned into the second pair. The difference is a total derivative that depends on the 
boundary properties. The last two terms, involving derivatives of M, are also known as Lifshitz invariants and were 
introduced by Mostovoy in the present context These terms are linear in P and therefore spontaneously brake the 
symmetry and make a nonzero P. In other words, as noted by Mostovoy, we have an effective local electric field Eint 
here, proportional to the magnetization gradients. 

E,nt = zsMiS/.M) + Z4[(M.V)A/] (34) 

Which leads to a polarization P = XsEmt- The overall effect back on the magnetization (within a local response) 
can be summarized as a free energy for magnetization: 

1 

F^- -Xij|z3M(V.M) + Z4[(M.V)Af]|2 (35) 

2XM 2 

The spatial profile of the ground state is determined by the momentum dependent XM{q), the q for which it is a 
maximum. The ferromagnetic ground state would be unstable if the effect of the Lifshitz invariant were to move the 
ground state to finite q. Note however that the effective free energy for M, resulting from the Lifshitz invariant is 
quartic in M. In other words, a ferromagnet is stable as long as 6m > {XE{q)zY ■ Since xe is divergent at T^;, a 
ferromagnetic ground state would be unstable near that temperature. 

The thermodynamic fluctuations, which refer to the contributions of inhomogeneous terms as well as higher order 
interaction terms, lead to a qualitative change in all temperature dependent properties. There are two effects: the 
first is a downward shift in the critical temperature. The second is the movement of asymptotic mean field behavior, 
resulting in a different exponent. The temperature range of these effects can be encapsulated into the Ginzburg 
parameters. 
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The fluctuation effects arise from gradient terms, 7e(VP)^ and 7m(VM)^, in tfie free energy (Eq. [22]) . Near Teq the 
corresponding correlation length diverges. The subleading length scale derived from 7™ starts to play an important 
role in the temperature range T^/ < T < Teo, but the relative role of 7™ and 7e changes as one moves between Teo 
and Tm- Near Tj\/ there is an effective free energy similar to Eq. I35[ which includes Z3 and 24, however a Gaussian 
calculation notes the effect of these terms as quartic in M and therefore negligible. These terms are also irrelevant in 
the RG sense in three dimensions, {x n> scales as b^^'^ in d-dimensions. 

The fluctuation contribution to the specific heat was computed in the gaussian approximation, and shown in Fig. 

El 
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FIG. 5: The most singular contribution to the specific heat from gaussian fluctuations, for different sets of correlation lengths 

We first include gradients, 7e(VP)^ and 7„i(VM)^ in the free energy. Then, we expand about the saddle point, 
mo+Sm, pQ+Sp, and neglect cross-terms since this is meant to be an estimate. There are two independent contributions 
in that approximation, one for the gaussian integral over the polarization and one from the magnetization. The most 
singular terms in the specific heat take the approximate form 

afiuct^T^G{T)F{T)^-^ 

Where the functions G and F follow a pattern. Let Oi denote the ith order parameter, M or P. Then G{T) = 
a/Tc + 120i^ + 2^02^ and F{T) = - 1) + 6biOl{T) + kOl{T). The correlation lengths will take the 
form, — -fj^ ■ The width of the critical region around either Tc in three dimensions-^ is given by the expression 

1/2 h 

TJ ~ 4AC jg? ' '^bere ACmf is the size of the jump in the specific heat from mean field theory. We plot the behavior 
of these decoupled fluctuations in Fig. [5l It will take a measurement to determine to what extent the system can be 
treated within mean field theory. We expect the dispersion of polar phonons to be relatively weak, and relying on 
experience in ferroelectrics, the polar transition could even be first order. 

It is instructive to consider the renormalization group flow and analyze the fixed points of our free energy. The 
critical behavior of two coupled order parameters each with 0{n) symmetry has already been computed in other 
contexts^^, and the results can be profitably reapplied in this new context. The necessary diagrams are shown in fig. 
ini The five fiow equations, resulting from integrating out a shell of momentum ^, A, are two copies of: 

'n = b^ln + 4(n, + 2)u, / + 2n,k / (36) 

b b 






FIG. 6: The diagrams at one loop, a) two-point function b) vertex for one order parameter c) vertex coupling two order 
parameters. Dotted or solid lines refer to the propagators for each order parameter separately. 



and the flow for the coupling of the two order parameters: 

'k = b'[k- 16k' / „ - 4(n,„ + 2)ku„, / ^ ^ - 4(n, + 2)ku, / ^ ^] (38) 

There are six fixed points determined by the fiow of the quartic terms Ui and k. In four of these, the order parameters 
decouple, k = 0. The remaining two non-trivial ones, the Heisenberg-Heisenberg and so-called biconical fixed points, 
interchange their stability as a function of how many components the order parameters have. Fluctuation driven 
first order transitions^^ can occur if there is a runaway flow in the RG (when there is no fixed point accessible under 
RG iterations), and are absent in this model. At the double- Heisenberg fixed point, Umc = Uec ^ kc ~ ^'^^n'^+s • 
Linearizing, the eigenvalues are all negative for Ug + rim < 4 and e > 0, but changes its stability if the components 
of the order parameter change. For the case n=3, that is for three dimensional vectors P and M considered here, 
the biconical fixed point is stable, and the double Heisenberg fixed point is unstable. If we have two easy planes 
instead and n=2, or very strong crystal field anisotropic such that n=l effectively, then the double-Heisenberg fixed 
point is again stable. The anomalous dimension as in all one-loop scalar field theories vanishes, but we can anticipate 
corrections at higher order. 



V. SUMMARY AND DISCUSSIONS 



We have studied the thermodynamic properties of a system with two order parameters and two transitions. Since, 
as a specific example, we have in mind a ferromagnet and a ferroelectric, the homogeneous interaction term is 
simply M'^P'. All the lowest order symmetry allowed interactions (without reference to crystallographic groups) 
were discussed in section IV. We calculate the phase diagram, order parameter temperature and field dependence, 
specific heat and susceptibility for the coupled order parameters. We follow through with a brief discussion of a 
renormalization group analysis. 
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Our results are (1) the cross susceptibility exists only below Tm, with a characteristic temperature dependence 
that diverges at Tm in the thermodynamic hmit. (2) The electric susceptibihty changes discontimiously at Tm- The 
temperature Tm{B,E) as a function of E follows a second order phase boundary. The effect of fluctuations (measured 
by the coefficients 7^ and 7e) would be to move T^o ^^nd Tm lower by amounts depending on 7^ and 7e. For large 
7m and 7,. the fluctuations play insignificant roles. The fiiictuation contributions in the vicinity of Tfjo are expected 
to be large. We expect though that qualitative features of this paper remain intact. The question remains whether 
the cross susceptibility is generated by fiiictuations at least in the temperature range Tm < T < Teq. There is no 
homogeneous factor that we can imagine which might lead to a nonzero Xi2- The possibility of a defect generated 
X12 remains a topic of future study. 
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